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Abstract 
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0^ ! The generalized Fermi-Dirac functions and their derivatives are important in eval- 

' uating the thermodynamic quantities of partiahy degenerate electrons in hot dense 

■ stellar plasmas. New recursion relations of the generalized Fermi-Dirac functions 

', have been found. An effective numerical method to evaluate the derivatives of the 

generalized Fermi-Dirac functions up to third order with respect to both degeneracy 
and temperature is then proposed, following Aparicio [14]. A Fortran program based 
on this method, together with a sample test case, is provided. Accuracy and domain 
of reliability of some other, popularly used analytic approximations of the gener- 



2 ' alized Fermi-Dirac functions for extreme conditions are investigated and compared 

c/3 . with our results. 
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Computers: Sun E4500/E5500, Compaq DEC Alpha, SGI Origin2000, HP 
Convex Exemplar, Cray SVl-lA/16-8, AMD K6 PC, IBM SP2 

Operating systems under which the program has been tested: Solaris 5.6, Linux 
(Red Hat 5.2), IRIX 64, SPP-UX 5.3, Unicos 10.0.0.6, Microsoft Windows 98 
(2nd ed.), AIX 4.2 

Programming language: FORTRAN 77 

Number of bytes in distributed program, including test data, etc.: 20758 

Distribution format: uuencoded (gzip) compressed tar file 

Keywords: Fermi-Dirac functions, equation of state: electron gas, numeric 
method, astrophysics plasma, stellar evolution 

Nature of physical problem 

Provide numerical method to evaluate generalized Fermi-Dirac functions and 
their derivatives with respect to r/ and f5 up to third order. The results are 
important for a highly accurate calculation of thermodynamic quantities of 
an electron gas with partial degeneracy and relatively high temperatures with 
very high order of accuracy. 

Method of solution 

Following the scheme proposed by Aparicio [|I^ , the generalized Fermi - Dirac 
integration is split into four optimized regions. Gauss-Legendre quadrature 
is used in the first three pieces, and Gauss-Laguerre quadrature in the last 
part when the e~^ term in the integrand dominates. Different break points are 
individually chosen for each t] derivative. 

Typical running time 

Less than 1 ms for each data point on a DEC Alpha station with a 533 MHz 
CPU in double precision. 

LONG WRITE-UP 



1 Introduction 

Relativistic quantum statistics is important for studying thermodynamic prop- 
erties of electrons at the elevated degeneracy and relatively high temperature 
as found in the interior of massive stars. Chandrasekhar expresses the 
thermodynamic quantities as integrals of hyperbolic functions [see Eq. (170) - 
(172) of |jl|]. In a more general way Cox & Giuli [^, write them in term of the 
generalized Fermi-Dirac (FD) functions with different orders [see Eq. (24.98) 
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- (24.100) of 0]. They define the generahzed FD functions as 





where t] = ^/ksT and (3 = k-BT/meC^ are the degeneracy parameter and a 
dimensionless temperature, respectively. Here, rUe is electron mass, T absolute 
temperature, /i the chemical potential, c speed of light in vacuum and k-Q the 
Boltzmann constant. The non-relativistic form of Eq. (1) (corresponding to 
/? = 0) is 

oo u 

i^T^^dx (2) 





is usually called the Fermi-Dirac function, and is used to study thermody- 
namic properties of the degenerate non-relativistic electron gas. For the usual 
physical applications, both for the non-relativistic and relativistic functions, 
the index k is an integer or half integer with k> —1. 

No exact analytic formulae are available to evaluate the generalized FD func- 
tions in full range of degeneracy and temperature. Some approximate expres- 
sions exist for some extreme cases, such as when degeneracy is very low, or 
extremely high 0; ^; ^ |^. There are also asymptotic expressions for the ther- 
modynamic quantities alone 0. On the other hand, Eggleton et al. 0, give 
a very useful global fitting formula for the generalized FD functions [see also 
[]5|], even though its accuracy is relatively modest. 



Numerical evaluation of the generalized FD functions is not an easy task. 
Simple tabular interpolations can cover only a small range of degeneracy 
and temperature, and the accuracy is not very good. Because the integrand 
is proportional to x'^ when x — 0, and decays as e~^ when x ^ oo, few 
simple numeric integration methods can simultaneously cover a sufficiently 
large range of degeneracy and temperature with acceptable accuracy as well 
as with an acceptable calculation effort. |jlO| uses a modified Gauss quadrature 



with weight function w{x) = x^/ (e^' +1), and [|Tl|] uses weight function w{x) = 
x^ I e^ to integrate the generalized FD functions directly. Both methods are 
relatively fast, but are only good for < 5. In |^, a Gauss-Laguerre quadrature 
with adjusted decay factor is used to cover — 1 < 77 < 11, and approximation 



formulae beyond that region. In [HI; O], different variations of the trapezoidal 



rule are being used to do the integration, but both involve a heavy trade-off 
between accuracy and computing time. 
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Mohankumar & Natarajan use another approximation method by sub- 
tracting the sum of residues of finite number of poles within a pre-selected 
contour from results of trapezoidal integration. This method is fast when rj is 
not too big, and the results can reach 14-digit accuracy in double precision for 
f3 < 100, and when rj is not too small, typically from rj > —25 to rj = 2000. It 
is good to about 10 digits when rj is as large as 10000. However, the number 
of poles to be included increases rapidly when rj gets too big. Extending this 
method to calculate different derivatives of the generalized FD functions with 
the same order of accuracy is not an easy task. 



Aparicio JT^ proposed another method which consists of splitting the generalized- 



FD-function integral into four regions. The break points are chosen such that 
errors from all four regions are balanced, and the location of the break points 
is optimized. Based on the property of the integrand, the Gauss-Legendre 
quadratures used are from x = to x = S3, and the Gauss-Laguerre quadra- 
ture from X = S3 to X — s> 00 [see Appendix (B)]. The method is fast (pro- 
portional to the number of points used in Gauss quadratures) and accurate. 
Extending it to include derivatives of the generalized FD functions is straight- 
forward. Our numerical method to evaluate derivatives of the generalized FD 
functions follows this scheme closely. We discuss it in more detail in Sect. 3. 

Besides the generalized FD functions themselves, their derivatives are also 
important for computing the thermodynamic quantities of an electron gas. 
Miralles & Van Riper |^ give the results of the first order of derivatives of 
the generalized FD functions with respect to both rj and (3. They use approx- 
imation formulae when t] is either very large or very small, and interpolate 
through tables when 77 is in between. The disadvantage is that it is only ac- 
curate to about 6 digits, and discontinuity exists when two different methods 
join. Also, their results are limited to k = |, | and | only. On the other hand, 
solar equation of state studies require up to third order derivatives on the 
generalized FD functions [0]. In particular, the precise data from helioseis- 
mology have demonstrated that the relativistic electron gas is necessary for 
solar modeling 0. For solar applications, high-precision computations of the 
relativistic thermodynamic quantities are mandatory . 



Because of its complexity, the generalized FD functions have not been fully 
understood. In Sect. 2 we first present some newly found recursion relations, 
and then discuss some of their properties. The numerical method to calculate 
derivatives of the generalized FD functions up to the third order is introduced 
in Sect. 3, based on the scheme of \Ti\. At the end we check the efficiency and 



validity of some popularly used analytic approximation formulae in several 
extreme cases in Sect. 4. 



4 



2 Recursion Relations 



Cox & Giuli show that the generahzed FD functions obey the following 
recursion relations 

P = {k + l)Fkir],P) (3) 



which link the FD functions Fk{ri,P) and Ffc+i(?7,/5) via both their r] and /? 
derivatives. Beyond that, we have found that the generalized FD functions 
also obey another recursion relation 

(4, 



which is between the generalized FD functions and their j3 derivatives. 

From Eq. (3) and Eq. (4) we can further derive two useful recursion relations, 
in which 



drj \ ' 2J ' dp 



is a relation among Fk-i{ri,P), Fk{ri,P) and Fk+i{ri,P), with both r] and P 
derivatives, and 

1-17 ( ^\^(^ r\ 9Fkiv,P) , 1^ dFk+iiv,P) 



is between the generalized FD functions and their t] derivatives only. 

These recursion relations can be used to check the accuracy of numerical 
evaluations as well as to simplify some theoretical expressions. However, one 
should keep in mind that when rj is very small the two terms on the right-hand 
side of either Eq. (3) and Eq. (5) are so close to each other that almost perfect 
cancellation happens. For this reason, the use of recursion relations to evaluate 
the derivatives of the generalized FD functions is not always possible. 

One interesting property of the recursion relation in Eq. (3) is that, for the 
non-relativistic limit case [P = 0), it can be written as: 

^^ = («^ + l)n.W (7) 
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Fig. 1. PLEASE PLACE FIGURE (1) HERE, 
which can be further written as: 



f X ^ f X X , 

/ QX / QX 

J e^-''+l J e^-'^+l (e^-^ + 1) 

\ /\ / 



Both Eqs. (7) and (8) are satisfied for arbitrary rj. This is very remarkable, 
because although the function x/[(e^~'' + 1)] is not equal to 1, when 
multiplied with the Fermi-Dirac kernel a;*'/(e^~'' +1), it leads to exactly the 
same value of the definitive integrals as if it were equal to 1. 

Another interesting behavior we observed from the generalized FD functions 
is the dependence on rj of the following ratio (see Fig. 1). 

{[dFu{r^,^)/dri]+P[dF,+,{^,P)/dri]Y 



The expression appears in the derivatives of number density or internal energy 



of the relativistic electron gas (see |jT6[ for more detail) 



For small degeneracy {rj — >■ — oo) the ratio (9) is very close to 1 for all value of 
index k {k > —1) and all (3 > 0. However, when the electron gas is extremely 
degenerate {t] oo), the ratio (9) approaches different limits for the non- 
relativistic and relativistic cases. For the non- relativistic case {(3 = 0), the 
limits of (9) are f ? |) |) ••• when /c = — respectively. The general 
expression is for k > —1 and k ^ 0. For the relativistic case (/3 7^ 0) the 
limits of (9) are |, |, |, ••• when A;=— respectively. The general 
expression is {k > —1). This result is true for 77 /5 1. So for a given 

small /3 (/? 7^ 0), that is, in the nearly non- relativistic case, for sufficiently 
large r], the relativistic relation is satisfied, as shown in Fig. (1). This property 
links the behavior of the generalized FD functions of both extreme cases (very 
small degeneracy or extremely great degeneracy). In Appendix (D) we will 
give a proof of this property using asymptotic expansions. 

One thing we want to point out is that, because the ratio of Eq. (9) is very close 
to unity when rj is small, one has to be cautious when evaluating the difference 
between the nominator and denominator for the case of a slightly degenerate 
electron gas in order to avoid numerical noise due to near cancellation. 
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Fig. 2. PLEASE PLACE FIGURE (2) HERE. 
3 Numerical Method and Test Run Results 



We choose the scheme proposed by [jT4| to calculate derivatives of the gener- 



alized FD functions, because this method is fast, accurate and transparent. 
The expressions for the derivatives are listed in Appendix (A). In the method 
adopted, the integration is further split into four domains. Gauss-Legendre 
quadratures are applied to the first three pieces, starting from x = to a pre- 
selected point S3, as defined in Appendix (B). In the last part, where x is from 
S3 to 00 the e~^ term dominates, and hence Gauss-Laguerre quadrature is 
used. The positions of the break points are optimized such that the combined 
errors from all the pieces contributing to the generalized FD functions are 
minimized and equally distributed among the pieces. Fortunately, the choice 
of break points turns out rather uncritical for each individual derivative [see 



also [14 



As shown in Appendix (A), the integrand of each derivative can be represented 



as a product of the generalized FD kernel x^^Jl + {l3x/2) /(e^~''+l) with a 
tuning function. When 77 and (3 are not very large, this tuning function varies 
mildly [see Fig. (2)]. Hence the same break points chosen for the generalized 
FD functions themselves can also be applied to their derivatives, and good 
accuracy is achieved. However, when 77 and (3 are getting large, new break 
points must be chosen. The 13 derivatives are not very sensitive to the choice 
of break points. Therefore, we provide new break points only for 77 derivatives. 
Formulae to calculate break points can be found in Appendix (B). 



One concern in the calculation of the generalized FD functions and their 
derivatives is how one determines the accuracy of the results. For the gen- 
eralized FD functions we compare the results of (the method adopted in 



the present paper) with those of |13[ evaluated in double precision. We find 
that they agree with each other to 14-digit accuracy for 77 up to 1000, and to 
10-digit accuracy for 7/ up to 10000. Additionally, both of them agree with [^ 
to 8-to-lO-digit accuracy for 77 up to 1000. 



For the derivatives of the generalized FD functions, we first test them using 
the different recursion relations available. Except for the aforementioned can- 
cellations when 1] is small, the recursion relations are satisfied with better than 
10-digit accuracy, for all derivatives with j3 being between 10~^ and 10'^ and 7/ 
between -100 and 10000. The only exception is the third-order derivative with 
respect to rj. Here, the same accuracy can only be achieved for rj less than 500. 
All of the above results refer to double precision (128 bits) (which is single 
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precision on a Cray). Of course, better results can be achieved with quadruple 
precision. In practice, if /c 7^ — |, the recursion relations can be used to convert 
the third-order rj derivative into a lower-order one when 77 is big, and in this 
way, very accurate results can be achieved. 

As a second test, we compare our directly evaluated derivatives with those ob- 
tained from numerical differentiation. For example, we compare dFi/2 {Vy P)/ 9r] 
from our method with Fi/2{rj + /?) — -^1/2(^7 — A?7, [3) /2 with the dif- 
ference Ar] chosen to be 10~^, 10~^, 10~®, 10~*, respectively, and the calcu- 
lation is done in double precision. The purpose of the test is (i) to test the 
accuracy of our method, and (ii) to find the dependence of Ar] as a function 
of accuracy. 

We conclude that both methods agree well. In practice, a step of Ai] or A/5 
between 10~^ and 10"'^ is an optimal choice to compute the derivatives of 
the generalized FD functions by numerical differentiation. For example, the 
numerical evaluation of the second-order derivative with respect to rj from the 
first-order derivative yields an accuracy of at least 8 digits in double precision. 

However, we would like to stress that our direct evaluations are superior to the 
numerical differentiation used here only for testing purposes. The accuracy of 
numerical differentiation can drop to half of the machine precision [i.e. to 8 
to 9 digits, compared to the 16 digit accuracy of double precision (non-Cray 
machines) or single precision (Cray machines)] when the aforementioned 3- 
point numerical differentiation is applied once. Applied twice, the resulting 
accuracy can drop to 3 to 5 digits for r] < 2000, and only 2 digits when r] 
5000. Applied three times, the resulting accuracy would only be 2 digits when 
t] < 50, and for higher rj values, the results would be completely useless. This 
demonstrates the well-known fact that numerical differentiation is a dangerous 
procedure if it is used beyond first-order derivatives. In addition, the step size 
to be used has to be chosen judiciously. For high-precision applications such 
as found in helioseismology, which requires precise second- and third-order 
thermodynamic quantities, the direct approach method presented in this paper 
is by far superior to more-than-one-time numerical differentiations. 

In this paper we have always used 200-point Gauss quadratures. Although 
fewer-points quadratures are also good for the calculations of the general- 
ized FD functions |T^, the 200-point Gauss quadratures used here guarantee 



the accuracy of the high-order derivatives. If computing speed is important, 
and third-order 77-derivatives are not needed, then the 40-to-80-point Gauss 
quadratures in WM are a good choice. 



To use the program accompanying this paper, the user has first to compile and 
link the three Fortran files. Then, upon choosing the index of the requested 
generalized FD function k, the label IB denoting either the function itself or 
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Table 1 

PLEASE PLACE TABLE (1) HERE. 



Fig. 3. PLEASE PLACE FIGURE (3) HERE. 



Fig. 4. PLEASE PLACE FIGURE (4) HERE. 



a specified derivative (see Table 1), the range of degeneracy t] and dimen- 
sionless temperature /?, the program then calculates the desired generalized 
FD function or derivative. A sample input is given in the file input, and the 
corresponding results are contained in the file testResult. 



4 Validity of Analytic Approximations 

Although the present numerical method to evaluate generalized FD functions 
is accurate within a very large range of input parameters rj and (3, analytic 
approximations still have advantages for theoretical studies and the under- 
standing of the underlying physics. However, since analytic approximations 
are usually only valid for restricted parameter regions, we must know them 
for each of them. In the following, we limit ourselves to the properties of 
the k = \ functions only if its behavior is representative for the other cases. 
Otherwise, we discuss the cases of index k = \ and ^ = §• 

First we check the accuracy of the non-relativistic (NR) FD function as an 
approximation to the generalized FD functions in the weakly non-relativistic 
case (/3 -C 1) 



From Fig. 3, we can see that when both r] and (3 are small the approximation 
is good to 6 digits, and the accuracy drops when either 77 or is getting big. 
Considerably better agreement (twice as many accurate digits, as shown in 
Fig. 4) can be reached by including one more term is 0: 



(10) 



Ffc(r/,/3)~Ffc(r7) + -/3Ffc+i(r/) 



(11) 



For the ultra- relativistic case (/? ^ 1) 0, results of the following approxima- 
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Fig. 5. PLEASE PLACE FIGURE (5) HERE. 



Fig. 6. PLEASE PLACE FIGURE (6) HERE. 



Fig. 7. PLEASE PLACE FIGURE (7) HERE. 



Fig. 8. PLEASE PLACE FIGURE (8) HERE, 
tion is shown in Fig. 5: 



'12) 



The figure shows that more than 6-digit accuracy can be obtained with this 
formula if /3 > 0, and rj is not too small, typically 77 > 0, and more than 
10-digit accuracy when (3 > 2. 

When ?7 > we can see that the asymptotic expansions by Pichon [see 
Appendix C.l, Eq. (Cl-Cll); also Cox & Giuh [§]] are very good for Fi(r],P) 
(Fig. 6), with better than 12-digit accuracy when 77 > 120. However, the 
results are only moderate for Fi{ri, (3) (Fig. 7) and F^ij], (3), in which case the 
accuracy is only better than 6 digits when 77 > 5000, and much less accurate 
when 7] is smaller. For convenience. We list the formulae for large and small rj 
in Appendix (C). 



When ?7 is small, the expression of the generalized FD functions in terms of 
modified Bessel functions [see Appendix C.2, Eq. (C12-C17); Cox & Giuli [§]] 
achieve very high accuracy. We have used Mathematica 4.0 (for Solaris) to 
evaluate the modified Bessel functions in double precision. Since the leading 
terms in Eq. (C.15) - (C.17) are comparable when f3 is small, cancellation 
effect would appear. A power series expansion expression [^ is therefore more 
appropriate. We have combined the methods and show their results in Fig. 8. 

We have also checked the quality of the fitting formulae by |§| [see Appendix 
C.3, Eq. (C21-C26)] (Fig. 9). Specifically, we have adopted the expression 
given by |Q (his Eq. 5a - 5c). These fitting formulae cover a much larger range 
of T] and P than those asymptotic formulae above (but with lower accuracy). 
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Fig. 9. PLEASE PLACE FIGURE (9) HERE. 



and they are adequate if accuracy is not critical, and temperature not too 
high. For instance, with these fitting formulae we can obtain 4-digit accuracy 
when 7] < 1000 and (3 < 10~^, and 2-digit accuracy when t] < 10000 and 
/3<3x 10~^. 



5 Conclusions 



In this paper we have investigated properties of the generalized FD functions 
in some detail, and we have found some useful new recursion relations. We pro- 
pose an accurate and relatively fast numerical method to evaluate the deriva- 
tives of the generalized FD functions, following the algorithm introduced by 



Tj] . Limitations of direct numerical differentiation are also discussed. Finally, 
we have obtained the range of validity of some popular analytic approxima- 
tions for the generalized FD functions. 
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A APPENDIX: Derivatives of the Generalized FD Functions 



Derivatives of the generalized FD functions up to third order with respect to 
7] and (3. 



dr] J e^-'?+l 1+e^-^ 







dp J e^-'?+l A + 2l3x 
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drf 
drjdp 
a/32 

d'Fk{v,(3) 

d'Fk{v,P) 

d'Fk{v,P) 
a/33 



f3x 



1 -e'?- 



1 



dx 



qx-v +1 4 + 2/3x 1 + ev-^ 



dx 



^x'Jl + f x' 



QX-ri^l (4 + 2/3x) 



dx 



X 



1 + if (1 _ e'?-^)2 - 2e^-^ 



(1 +e^- 



x\3 



X 



e^-'?+l (1 + e''-^)2 4 + 2/3x 

1 



dx 



dx 



xS/l + f x2 



e^-'7+l (4 + 2/3x)2 1 + 



dx 



°^xVl + f 3x3 
e^-''+l (4 + 2/3x)- 



■ dx 



B APPENDIX: Break Points 



Three break points 5*1, 5*2, 5*3 are [as shown in |M[ ]: 



•S*! = Xa — Xb 

52 = Xa 

53 = Xa + Xc 



where 



^ ai + 61^ + ci^2 

_ 02 + + C2<j2^^ 

l + e2e + C2e 
1 + ese + cs^ 



12 



^(r^) = (T-Mn[l + e'^^^-^)] (B.7) 
where the corresponding parameters are: 



Parameter 




drj 




drp 


D 


3.3609E0 


4.99551E0 


3.93830E0 


A 


1 7zLzLzLF,n 


a 


9.1186E-2 


9.11856E-2 


9.11856E-2 


Q 


± ± O 5j \ J J-l; 




6.7774E0 


6.77740E0 


6.77740E0 


u 


7774nEn 


hi 


1.1418E0 


1.14180E0 


1.14180E0 


1 

1 




Ci 


2.9826E0 


2.98255E0 


2.98255E0 
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02 


3.7601E0 


3.76010E0 


3.76010E0 


o 


7801 DRD 


b2 


9.3719E-2 


9.37188E-2 


9.37188E-2 


Q 


O 1 XOO-L; Zi 


C2 


2.1064E-2 


2.10635E-2 


2.10635E-2 


2 


10635E-2 


d2 


3.1084E+1 


3.95015E1 


3.14499E1 


3 


05412E1 


62 


1.0056E0 


1.00557E0 


1.00557E0 


1 


00557E0 


03 


7.5669E0 


7.56690E0 


7.56690E0 


7 


56690E0 


^3 


1.1695E0 


1.16953E0 


1.16953E0 


1 


16953E0 


C3 


7.5416E-1 


7.54162E0 


7.54162E0 


7 


54162E0 


d3 


6.6559E0 


7.64734E0 


6.86346E0 


7 


88030E0 


ea 


-1.2819E0 


-1.28190E-1 


-1.28190E-1 


-1.28190E-1 



C APPENDIX: Asymptotic Expressions for Great and Small De- 
generacy 

C.l Great degeneracy (rj > 0) /0; H/ 



Fi(r/,/3) = -= {f,^{y) + {l + r^p)[C, + C2 + C,{Ay' + 7)]} (C.l) 



F| iv, /?) = {/| (y) + Ci (3 + 2.^/3) - 



-Cs[(Arj(3 + 6)rj(3 + 3]} 



(C.2) 



Fa (77, f3) = (/5 (y) + Ci(5 + 77/3) + C2{[(277/3 + 10)rjP + 15]vP + 5} 

+6-3(3 + 577/?)) (C.3) 



where 



7r2 



7 ^ /tt/?' " 



Ci — — , C2 — —Ci — 



6 1/ 



20 ^Vr 



31 



168 V 2/ 



(C.4) 
(C.5) 



For ri > 0.05: 



yyl + y^ — arcsinh(y) 



= (1 + ^y^) + - - ^arcsinh(y) 



(C.6) 
(C.7) 
(C.8) 



For rj < 0.05: 



{[( 



58773 



2 , 77 \ 2 "-^ 1 2 I 1 2_ 
^ ."^i^n/ ^ .^'^98J ^ 1408-1 ^ 



63 



1114112 " 5120> 



3328 



,2 I 35 -1 2 ^ \ „.2 



1408- 



144. 



5 ■ 
+36. 



(C.9) 



{[( 



58773 
1114112 



77 
5120 



1/^ + 



63 



3328 



35 ~| 9 5 



1408 



144. 



5 

36 



V + — 
^10 



(C.IO) 



275913 
4456448 



^ 512^ ^ 



416 



15 



704 J 



1 

36 



+ 



28j 



(C.ll) 
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C.2 Small degeneracy (^7 < Oj /j^; |^ 



Define K^{x) = e^ Ky{x) where Ky is the modified Bessel function, then: 



( j^^n— 1 



E 



1 °o ('_i\n-l C 



HI -HI 



(C.13) 

.fA,(;i)}(c.u, 



While for rj > —30 and log/5 < —1, the following series expansions by Miralles 
& Van Riper M are better: 



F,(r^,/5) = vlZ^e" fl + -/3 - — /^^ + i^^^/^^ - 
2^"^' 2 V 8 128 1024^ 1024^ 



vr 



Fs (r/, /?) = ^ e" 1 + -(3 
2^''^^ 2 V 8'^ 4090 



8' 128' 16384' 
539 ^, 



(C.15) 
(C.16) 
(C.17) 



C. 3 Global fit from EFF |/ 



If one changes the parameters of the generahzed FD functions from (r^,/?) to 
(/, 5() where: 



7] = 2m + In 



u- 1 



u+l 



P = g/u 



(C.18) 
(C.19) 



and 



(C.20) 



then the generahzed FD functions can be approximated as 



(C.21) 
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(C.22) 
(C.23) 



where 



P 



I _|_ \?,/2 l^m=0 l^n=0 PmnJ U 



(C.24) 



P = 



I q^/^n _)_ (7)3/2 Z^m=0 ^n=0 -^mnj U 



(C.25) 



M = 



i q^/^fl + q)^/^ Z^m=0 Z^ra=0 ^mnj 



(C.26) 



The values of the coefficients pmn, Pmn and f/mn can be found in p| with 
M = = 1, 2, 3 or 4. 



D APPENDIX: Proof of Eq. (9) 

We use the series expansion approximation to prove the following expression, 
for which k > —1. 



[F,{rj,P)+PFk+iiv,P)] {[d^F,ir^,/3)/dv']+(3[d'Fk+,iv,P)/dv^]} 

1 r] ^ — oo 
^ Mi ^_^^^^ = o,A;^0 (D.l) 



When 7] — oo, (3 = 0, the approximation of the generalized FD functions is 
given by Cox & Giuli || Eq. (24.111) as: 




oo,/3 ^ 



D.l Small degeneracy, non-relativistic 



Ffc(-cx), 0) = Ffc(-oo) = e" T{k + 1) 



(D.2) 
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So we get: 



Hence Eq. (D.l) can be written as: 

[dFu{-^)/drf 



(D.4) 



Fk{-^WFu{-^)ldr]^] 
D.2 Small degeneracy, arbitrary relativistic 



When rj ^ —oo, (3 ^ 0, we can find the leading term of the generalized FD 
functions for the following two extreme cases from ^ Eq. (24.253) - (24.255) 
and Eq. (24.267) - (24.269), respectively, as: 



Fk{-oo,P)=e^gk{P) /3 < 1 
F,,(-oo,/3) =e^/ife(/3) /3>1 



(D.5) 



For which we can see that: 



dF,{-^,(3) d'F,{-^,(3) 
F.i-oo, P) = = (D.6) 



Hence Eq. (D.l) is satisfied. 

D.3 Great degeneracy, n on- relativistic 

When r/ ^ cx), /3 = we get from g] Eq. (24.116) that 

F,(oo, 0) = Ffc(oo) = v'-' (D-7) 



So 



[dFkioo)/drif k + 1 



F,(oo)[92Ffc(oo)/9r/2 



(D.8) 
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D.4 Great degeneracy, arbitrary relativistic 



When rj oo, /? 7^ 0, the leading term of the generahzed FD functions from 
B Eq. (24.179) - (24.181) is: 



where = {l+r](3)'^. When r/ ^ cx), /5 7^ 0, y r^/? > 1. So, y/1 + y'^-1 
y. From || Eq. (24.202), we can also find 

So, 



Hence we can see that 

[dF,,{oo,(3)/d7]]^ fc + 3/2 

Fk{(^,(3)[d^Fk{(^,f3)/dr]^]^ k + 1/2 



(D.12) 



If we look at Eq. (D.l) we can see that when 77 00 the Fk+i{oo,(3) terms 

dominate over the Fk{oo,(3) terms. As a result, the ratio in Eq. (D.l) will be 

fc+5/2 

fe+3/2' 
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Table 1: Index of the function or derivative to be evaluated by the program. 



TR 
±±j 


J- LliU^. 


TR 

J.J_J 


J. LlliL^. 


n 

U 


^k\JI-i P) 


1 
± 


dFk(v,0) 
8r) 


2 


80 


3 


d^Fk{n,0) 


4 




5 


d^Fk{v,0) 
80^ 


6 


d^Fk{v,0) 


7 


d^Fkiv,0) 
d'rpd0 


8 


d^Fk{n,0) 


9 


d''F„{'n,0) 
803 
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Fig. 1: Different limits approached in Eq. (9). a,, k — h. k — ^, c. k — 
d. k — |. Solid line for /3 = 1, and dashed line ior P — 0. 



3 
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Fig. 2: Behavior of the tuning functions for (3 — 1 and 77 = 10 
case. Indexes refer to the corresponding derivatives in Table (1). 
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Fig. 3: Effective range of non-relativistic approximation for 
Fi{r],(3) from Eq. (10). y axis is the logarithm of relative differ- 
ence log{\Fapprox — F\/F); the same notation for Figs. (4) to (9). 




logp 



Fig. 4: Effective range of non-relativistic approximation 
Fi{r),P) from Eq. (11). Notations are the same as in Fig. 
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Fig. 6: Effective range of extremly-degenerate approximation for 
Fi{r],(5) from Eq. (C.l). Notations are the same as in Fig. (3). 




Fig. 7: Effective range of extremly-degenerate approximation for 
F3(ri,/3) from Eq. (C.2). Notations are the same as in Fig. (3). 



X 10^ 
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Fig. 8: Effective range of non-degenerate approximation for 
Fi{r),P) from Eq. (C.12). Notations are the same as in Fig. (3). 
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Fig. 9: Effective range fitting formula of Fi{ri,(3) by 
(see 1^). Notations are tlie same as in Fig. (3). 
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